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Recent work has shown that scalar fields around black holes can form long-lived, quasistationary 
configurations surviving for cosmological timescales. With this requirement, scalar fields cannot be 
discarded as viable candidates for dark matter halo models in galaxies around central supermassive 
black holes (SMBH). One hypothesis for the formation of most SMBHs at high redshift is the 
gravitational collapse of supermassive stars (SMS) with masses of ~ 10 5 Mq. Therefore, a constraint 
for the existence of quasi-bound states of scalar fields is their survival to such dynamic events. To 
answer this question we present in this paper the results of a series of numerical relativity simulations 
of gravitationally collapsing, spherically symmetric stars surrounded by self-gravitating scalar fields. 

We use an ideal fluid equation of state with adiabatic index T = 4/3 which is adequate to simulate 
radiation-dominated isentropic SMSs. Our results confirm the existence of oscillating, long-lived, 
self-gravitating scalar field configurations around non-rotating black holes after the collapse of the 
stars. 

PACS numbers: 95.30.Lz, 95.30.Sf, 04.70.Bw, 04.25.dg 


I. INTRODUCTION 

The formation of galactic black holes is a fundamental 
issue in the study of galaxy formation. Most galaxies 
host black holes in their nuclei. Very luminous quasars 
at redshifts higher than 6 are believed to be powered 
by supermassive black holes (SMBH) with masses of ~ 
1O 9 M 0 [2 |2J. However, it is still unknown how such 
SMBH can form and grow to reach the current values 
when the Universe was less than 10 9 years old. 

One possible scenario for the formation of SMBH in 
the center of galaxies is that they were formed through 
accretion and merging of the first stellar remnants E3HI- 
Another scenario assumes that central black holes formed 
via direct collapse (see EHZ]). This route assumes that 
fragmentation, which depends on efficient cooling, is sup¬ 
pressed. Depending on the rate and efficiency of the in¬ 
flowing mass, there may be different outcomes. A low 
rate of mass accumulation would favor the formation 
of isentropic SMSs, with masses > 5 x 1O 4 M 0 , which 
then would evolve as equilibrium configurations domi¬ 
nated by radiation pressure before collapsing to form a 
SMBH seed [8] 9]. 

Scalar fields have been considered in cosmology to pro¬ 
vide inflationary solutions in the early universe m and 
as an alternative explanation for dark energy HU- In ad¬ 
dition, they have also been proposed as strong candidates 
for dark matter [IMS]- In these models one assumes 
that dark matter is composed by bosonic particles which 
may condensate into macroscopic objects forming boson 
stars or forming halos around black holes. The fact that 
most galaxies host a SMBH at their centers has raised 


many questions regarding the dynamics of the whole sys¬ 
tem. 

The issue of whether scalar field halos are able to sur¬ 
vive around black holes for cosmological timescales has 
been addressed in several occasions (see e.g. [16] and ref¬ 
erences therein). It has been shown in the test-field limit 
in the background spacetime of a Schwarzschild black 
hole, that scalar field configurations do remain around 
SMBHs for cosmological timescales, provided the mass 
of the scalar particle is of the order of 10~ 22 — 10 _24 eV. 
Further studies considered also the backreaction of space- 
time dZHU], and confirmed the existence of such long- 
lived states. These gravitationally trapped, scalar field 
configurations are states characterized by precise com¬ 
plex frequencies. It has been proved that these states 
are generic in the sense that they appear after some 
time in the evolution of the scalar field even when ar¬ 
bitrary initial data is used [20]. This observation seems 
to be enough to guarantee (as we show below) that any 
choice of initial data will create a scalar held environment 
around the black hole. 

The main objective of this work is to assess the role of 
black hole formation in the presence of scalar held config¬ 
urations which resemble the quasi-bound states found in 
previous investigations where the black hole is assumed to 
exist from the start. To accomplish this goal we perform 
numerical relativity simulations which track the gravita¬ 
tional collapse of a polytropic star in the presence of a 
scalar held environment. The numerical techniques and 
the code used are those already described in [17] ■ We 
shall refer to that work as Paper 1. In this work we also 
take into account the evolution of the relativistic hydro- 
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dynamic equations. In order to simplify our model, we 
consider that the scalar held and the matter of the star 
are coupled only through gravity. 

The paper is organized as follows: In Section [TT] we 
present the basic equations and discuss the initial data 
used in our simulations. Section III briefly describes our 
numerical approach. In Section IV we discuss our find¬ 
ings and describe some properties of the solutions. Fi¬ 
nally in Section [V] we sum up our concluding remarks. 
Throughout the paper Greek indices run over spacetime 
indices (0 to 3), while Latin indices run over space indices 
only (1 to 3). We use units in which c = G = Mq = 1. 


II. BASIC EQUATIONS 

We investigate the dynamics of a self-gravitating scalar 
held conhguration around a polytropic star collapsing to 
a black hole by solving numerically the coupled Einstein- 
Klein-Gordon system and the general relativistic hydro¬ 
dynamics equations. The Einstein equations in covariant 
form read 

Rap ^ 9apR — $7rT a ^ ; (1) 

where R a p is the Ricci tensor of the 4-dimensional space- 
time, g a p is the spacetime metric, R is the Ricci scalar, 
and the matter content is given for our system by a com¬ 
bined stress-energy tensor 

T a p=Tlp+T%, (2) 

where T^g is the stress-energy tensor of the (perfect) huid 
dehned by 

T ap = phu a Uf3 + pg a p , (3) 

and is the part corresponding to the scalar held, 
given by 

= d a <f>dp$ - l -g a p + g 2 ® 2 ) . (4) 

In Eq. © p is the rest-mass density, p is the pressure, 
u a is the huid four velocity and h = 1 + e + p/p is the 
relativistic specihc enthalpy, with e being the specihc in¬ 
ternal energy. We will assume that each type of matter 
obeys their own conservation law, in other words, that 
the coupling between them comes purely from gravity. 
The equations of motion for the scalar held are obtained 
from the Klein-Gordon equation 

□ $ - g 2 <$> = 0 , ( 5 ) 

where the D’Alambertian operator is dehned by □ := 
(l/y/^g)d a (y/ zr gg a Pdp). We follow the convention that 
$ is dimensionless and g has dimensions of (length)^ 1 . 

In the following we present the different systems of 
evolution (and constraint) equations we solve to carry 
out our study. While we mainly include this information 
to make the paper self-contained, we keep these sections 
as concise as possible, and refer the interested reader to 
our previous papers HU [211121 for further details. 


A. Einstein’s equations 


Under the assumption of spherical symmetry the 1+1 
line element may be written as 


ds 2 = e 4x (a(t, r)dt 2 + r 2 b[t, r)dfl 2 ) , ( 6 ) 


with dtt 2 = sin 2 ddp 2 + dd 2 being the solid angle element 
and a(t, r) and b(t , r) two non-vanishing metric functions. 
Moreover \ is related to the conformal factor ijj as ij) = 
e x = ( 7 / 7 ) 1 / 12 , with 7 and 7 being the determinants of 
the physical and conformal 3-metrics, respectively. They 
are conformally related by 7 ij = e 4xx /ij. 

In particular we adopt the BSSN formalism of Ein¬ 
stein’s equations [231 121] , where the evolved fields are the 
conformally related 3-dimensional metric, the conformal 
exponent X) the trace of the extrinsic curvature K , the 
independent component of the traceless part of the con¬ 
formal extrinsic curvature, and the radial component of 
the conformal connection functions (we refer to [231 US] 
for further details). 

The matter fields appearing on the right-hand-side 
of the evolution equations for the gravitational fields 
(see [23123), £, ji, Si, and S,j include the contribution 
of both the fluid and the scalar field, i.e. £ = £ F + £ SF , 
etc. The matter source terms for the fluid given by the 
stress-energy tensor of Eq. © read 


ri a n p Tlp, 

(7) 

-7 > P Tlp, 

( 8 ) 

c T F y r , 

(9) 

C T F ) e e■ 

( 10 ) 


The Hamiltonian and momentum constraints are given 
by the following two equations that we compute to mon¬ 
itor the accuracy of the numerical evolutions: 

n = R~{A 2 a + 2A 2 b ) + ^K 2 -16n£ = 0, (11) 
2 

M r = d r A a - - d r K + 6 A a d r x 

+ (A a — A b )^-~ |—— 87T+-= 0, (12) 

where A a and A b are the contraction of the traceless 
part of the conformal extrinsic curvature, A a = A r r , A b = 
A e g . In addition to the evolution fields, there are two more 
variables left undetermined, the lapse function a, and 
the shift vector fi l . Our code can handle arbitrary gauge 
conditions and for the simulations reported in this paper 
we use the same conditions of Paper 1, namely the so- 
called “non-advective 1 +log” condition [26] for the lapse, 
and a variation of the “Gamma-driver” condition for the 

shift vector [23] [2Z1 • 
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B. Klein-Gordon equation 


We then assemble these variables into a vector U of con¬ 
served fluid variables 


To solve the Klein-Gordon equation we introduce two 
first-order variables defined as: 


U = ^7 (D,S r ,r). (27) 


n := n“<9 a $ = -(<9 t $-/3 r <9 r $) , (13) 

a 

'L := d r <S> . (14) 

Therefore, using Eq. © we obtain the following system 
of first-order equations: 

<9 t <J> = 
dtV = 
d t U = 

+ 

+ 

The matter source terms for the scalar field read 

s sp ^ + ( 18 ) 

if = -7= - 5 „.nvl/, (19) 

sf = ( T SF y r = l(n 2 + , ( 20 ) 

sf = (T SF )S = l(n 2 - . ( 21 ) 

C. Hydrodynamics equations 

The general relativistic hydrodynamics equations, ex¬ 
pressed through the conservation laws of the stress- 
energy tensor T and the continuity equation, are 


p r d r § + an , 
/3 r d r ^ + $d r /3 r 


9 r (an) 


S r d r u 


a 


ae 


4 X 


[5 r T 


'2 9 r a drb 


m 4 x 


9 r a + qAT — a/r 2 $ . 


(15) 

(16) 


(17) 


By defining the corresponding vector of fluxes F r and 
the vector of source terms S we can cast Eqs. (22) in 
flux-conservative form 


<9 t U + <9 r F r = S. (28) 

We refer to m where the explicit form of vectors F r and 
S is given. To close the system of hydrodynamics equa¬ 
tions, we choose a Gamma-law equation of state (EOS) 

P = (T-l)pe, (29) 


where T is the adiabatic index. 


D. Scalar field initial data 


As in Paper I our choice of initial data for the scalar 
field is a Gaussian distribution of the form 

$ = A 0 e " (7 - ro)2/A2 , (30) 

where A 0 is the initial amplitude of the pulse, r 0 is the 
center of the Gaussian, and A is its width. The auxiliary 
first-order quantities are initialized as follows 


n(t = 0,r) = 0 , (31) 

T(f = 0, r) = -2 (r ~ ro) A 0 e - (r - ro)2/A2 . (32) 

A 4 


Likewise, we choose a conformally flat metric with 
a = b = 1 together with a time symmetry condition 
Kij = 0. With these choices and II(t = 0, r) = 0 the mo¬ 
mentum constraint is satisfied trivially and the Hamilto¬ 
nian constraint, Eq. (11), yields the following equation 
for the conformal factor ip = e x , 


V a T af) = 0 , V« (pu a ) = 0. (22) 


Following [2S], we write the equations of general rela¬ 
tivistic hydrodynamics as a first-order, flux-conservative 
system. We define the fluid 3-velocity as seen by a normal 
(Eulerian) observer as 


au L 


f3 r 

a 


(23) 


and the Lorentz factor between the fluid and the normal 
observer as W = avP. We also define the relativistic 
densities of mass, momentum, and energy for the normal 
observer as 


D = pW , (24) 

Sr = phw 2 v r , (25) 

t = phW 2 -p-D. (26) 


d rr ip H— d r ip + 2mp 5 £ = 0 . (33) 

r 

In order to solve this equation for our collapsing stars, 
we assume that the conformal factor can be written as 

iP = ip T °v + u (r) , (34) 

where ip TOV is the part of the conformal factor as¬ 
sociated with the stellar Tolman-Oppenheimer-Volkoff 
(TOV) model. By substituting this ansatz in the Hamil¬ 
tonian constraint we obtain 

d rr u(r)+ 2 d r u(r)+2nip 5 £ SF +2n{ip 5 -(ip TOV ) 5 )£ TOY = 0 . 
r 

(35) 

Given a distribution of the scalar field density £ SF , we 
solve this ordinary differential equation for u using a 
standard fourth-order Runge-Kutta integrator, assuming 
that u —» 0 as r —► oo and regularity at the origin. 
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TABLE I: Stellar TOV models: N = 3 (T = 4/3) polytropes 
with k = 0.48773. From left to right the columns indicate the 
model number, the central rest-mass density, p c , the central 
energy density, e c , the radius of the star, R, the gravitational 
mass of the star, M, and the compactness, M/R. 


Model 

Pc 

e c 

R 

M 

M/R 

1 

3.11e-6 

3.17e-6 

185.8 

1.43 

0.008 

2 

3.11e-5 

3.25e-5 

86.3 

1.31 

0.015 

3 

3.11e-4 

3.41e-4 

39.9 

1.01 

0.025 

4 

3.11e-3 

3.77e-3 

19.2 

0.80 

0.042 



FIG. 1: Mass of the polytropic TOV stars as function of their 
initial central rest-mass density. All models are in the unsta¬ 
ble branch of the diagram. 


III. NUMERICS 

The time update of the different systems of evolution 
equations we have to solve in our code (Einstein, Klein- 
Gordon, and Euler) is done using the same type of tech¬ 
niques we have extensively used in our previous work 
(see in particular [TT] [5Tj). We address the interested 
reader to those references for full details on the particu¬ 
lar numerical techniques implemented in the code. Here, 
we simply mention that the evolution equations are inte¬ 
grated using the second-order PIRK method developed 
by [29] 130) . This method allows to handle the singu¬ 
lar terms that appear in the evolution equations due to 
our choice of curvilinear coordinates. The derivatives 
in the spacetime evolution are computed using a fourth- 
order centered finite difference approximation on a uni¬ 
form grid except for advection terms for which we adopt 
a fourth-order upwind scheme. We also use fourth-order 
Kreiss-Oliger dissipation to avoid high frequency noise 
appearing near the outer boundary. The difference with 
our previous work HZ! is that we are also evolving the 
hydrodynamics conserved quantities explicitly. Corre¬ 


spondingly, the equations of hydrodynamics are solved 
using the HLLE approximate Riemann solver in tandem 
with the second-order MC reconstruction scheme m- 


IV. RESULTS 

A. Stellar models and code convergence 

We construct spherically-symmetric stellar models by 
solving the Tolman-Oppenheimer-Volkoff (TOV) equa¬ 
tions for a polytropic EOS, 


p= K p 1+1 / N , 


(36) 


where k is the polytropic constant and N the polytropic 
index. These data are evolved in our code with the 


Gamma-law EOS, Eq. (29), with T = 1 + 1/iV. We adopt 
N = 3 (T = 4/3), which is the appropriate value for a 
SMS. Moreover, we select the value of n = 0.48773 in our 
code units. The polytropic constant is set to this value 
to keep the mass of the star close to one, namely in the 
range M = {0.8,1.43} for our four different TOV mod¬ 
els. These models are gravitationally unstable spherical 
relativistic stars that collapse to black holes in dynami¬ 
cal timescales by artificially decreasing the pressure by a 
small amount (typically 5%) at t = 0. The main charac¬ 
teristics of our sample of stellar models are reported in 
Table [T] Their corresponding location in the mass vs cen¬ 
tral density diagram is displayed in Fig. |Tj The models 
are constructed using the RNS code [51] . 

The collapsing polytropic stars are surrounded by a 
scalar field cloud whose functional form is given by 
Eq. (30). We investigate the behaviour of the scalar 


field in a highly dynamical evolution in which the stars 
undergo gravitational collapse. The end result of the 
collapse is the formation of a Schwarzscliild black hole 
whose mass is the sum of the mass of the corresponding 
original star and a portion of the scalar field mass, the 
part that is accreted. The aim of this investigation is to 
find out whether in such a dynamical scenario, part of 
the scalar held can still survive in the form of long-lived 
quasi-bound states. 

To carry out our study we set up and evolve 80 dif¬ 
ferent initial scalar held models. These models are char¬ 
acterized by different values of the initial energy of the 
scalar held, namely E 0 = {3.2 x 10 -5 ,0.50,2.85,4.3} 
and different values of the scalar held mass p = 
{0.05,0.08,0.1,0.15,0.2}. The choice of the values for 
the energy is made to span all regimes of interest, from 
the test-held regime to that for which the self-gravity of 
the scalar held becomes important. For all models, the 
distribution of the scalar held pulse is placed at the same 
initial location = 100 and they all share the same ini¬ 
tial width A = 50. 

We have used a logarithmic radial grid in our simula¬ 
tions setting to A r = 0.1 the minimum radial resolution 
at the origin. This is different from the grid we used in 
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FIG. 2: Extent of the radial grid as a function of the num¬ 
ber of grid points. The black solid line represents the radial 
coordinate profile of the composite grid, the red dotted line 
the geometric progression patch and the green dashed line the 
hyperbolic cosine patch. 



FIG. 3: Time evolution of the L2 norm of the Hamiltonian 
constraint for an initial scalar field pulse with Eq = 5 x ICG 3 
and mass Mfi = 0.1 for the TOV model 4. The plot shows 
results for three different resolutions, rescaled by the factors 
corresponding to second-order convergence. 


Paper 1. The reason behind this choice is to reduce the 
number of radial points without loosing accuracy in the 
inner region and placing the outer boundary sufficiently 
far away to avoid unwanted reflections from the scalar 
field. Our grid is composed of two patches, a geometrical 
progression in the interior and a hyperbolic cosine out¬ 
side. Using the inner grid alone requires too many grid 
points to place the outer boundary sufficiently far from 
the origin, while the hyperbolic cosine patch produces 
very small grid spacings in the inner region of the grid, 
leading to prohibitively small timesteps due to the CFL 
condition. The two patches are defined at each radial 
point as 


r'geo(^) — ^ llprog T’min (37) 

r c osh(«) = cosh[i (—log((l + ^max) 

+ \/2r max + ^max)/(7^r))] — 1) (38) 

where A?' is the minimum resolution, r m i n and r max the 
inner and outer boundary respectively, and N r the to¬ 
tal number of radial points. The outer boundary of the 
computational domain is placed at r max = 4x 10, far 
enough as to not affect the dynamics in the inner region 
during the entire simulation. For the matching of the two 
patches, we compute the grid spacings Ary for the two 
patches r geo (i) and r cos h(*) at each grid point and when 
the difference between the grid spacings is at a minimum, 
we change from the interior grid to the hyperbolic cosine 
grid from the corresponding radial point onwards. In 
fig- 0 we plot the two patches and the resulting grid as 
a function of the number of grid points. The time step is 
given by At = 0.5A?’ in order to obtain long-term stable 
simulations. The final time of the numerical evolutions 


is 3.5 x 10 4 . 

While the convergence of the code was already satis¬ 
factorily tested in Paper 1, we nevertheless need to assess 
it again here since we now deal with an augmented sys¬ 
tem of equations than in Paper 1 due to the presence of 
the hydrodynamical terms. In order to test the conver¬ 
gence of the code we performed three simulations with 
different resolutions Ar = {0.2,0.1,0.05}. In Fig. [3] we 
plot the rescaled evolution of the L2 norm of the Hamil¬ 
tonian constraint for a particular choice of the initial en¬ 
ergy, E 0 = 4.3, scalar field mass, /x = 0.1, and for the 
TOV model 4, obtaining the expected second-order con¬ 
vergence of our PIRK time-evolution scheme. We note 
that the same result is achieved irrespective of the combi¬ 
nation of parameters considered and of the TOV model. 

Comparing simulations with different resolutions, we 
estimate that the error in the final black hole mass is 
~ 2% while the error estimates for the scalar field en¬ 
ergy and ADM mass at the end of the simulations are 
~ 2 — 3%. Those errors decrease as (A?-) 2 with increas¬ 
ing resolution. Taking into account that our simulations 
are performed up to a (significantly long) final time of 
3.5 x 10 4 , the error bars for our main quantifiable quan¬ 
tity, the frequencies of oscillation of the scalar field (see 
next section), are ~ 3 x 10~ 5 . 


B. Non-linear quasi-bound states 


We solve the Einstein-Klein-Gordon-Euler system us¬ 
ing the initial data given by Eqs. (30)-(32) and let the 


scalar field evolve in the rapidly changing gravitational 
field of a collapsing polytropic star. As in Paper 1 we an- 


















TABLE II: Initial parameters and most relevant quantities for the scalar field and collapsing TOV model 1 (M = 1.43). From 
left to right the columns report: the name of the model, the initial amplitude of the pulse, Ao, the scalar field mass, /r, the real 
part of the angular frequency uj for the fundamental mode of oscillation and the first overtone, the time of the collapse of the 
TOV star and the formation of the black hole, t c 0 i, the final mass of the apparent horizon of the black hole, Mah, the initial 
ADM mass of the system, Madm, and the initial and final energy of the scalar field, Eo and Ea na i respectively. The initial 
Gaussian pulse is located at r o = 100 with half-width A = 50. 


Model 

^0 

A 1 

CO 

1 2 

^col 

Mah 

Madm 

E 0 

-S'final 

l_la 

5.00E-5 

0.05 

0.04990 


2494 

1.40 

1.43 

3.2E-5 

3.1E-5 

l_lb 

3.30E-5 

0.08 

0.07935 

0.07988 

2494 

1.40 

1.43 

3.2E-5 

2.0E-5 

l_lc 

2.65E-5 

0.10 

0.09874 

0.09964 

2494 

1.40 

1.43 

3.2E-5 

1.2E-5 

l_ld 

1.80E-5 

0.15 

0.14901 

0.14954 

2494 

1.40 

1.43 

3.2E-5 

1.2E-5 

l le 

1.35E-5 

0.20 

0.19783 

0.19927 

2494 

1.40 

1.43 

3.2E-5 

5.0E-6 

l_2a 

6.30E-3 

0.05 

0.05009 


1706 

1.40 

1.92 

0.50 

0.49 

l_2b 

4.12E-3 

0.08 

0.07953 

0.08007 

1582 

1.67 

1.92 

0.50 

0.22 

l_2c 

3.33E-3 

0.10 

0.09859 

0.09999 

1508 

1.76 

1.92 

0.50 

0.15 

l_2d 

2.25E-3 

0.15 

0.15008 


1376 

1.77 

1.92 

0.50 

0.14 

l 2e 

1.70E-3 

0.20 

0.20017 


1306 

1.90 

1.92 

0.50 

0.02 

l_3a 

1.50E-2 

0.05 

0.05099 


774 

1.49 

4.20 

2.85 

2.71 

l_3b 

9.82E-3 

0.08 

0.07909 

0.08115 

748 

3.38 

4.19 

2.84 

0.81 

l_3c 

7.95E-3 

0.10 

0.10197 


738 

3.47 

4.19 

2.84 

0.70 

l_3d 

5.36E-3 

0.15 

0.15314 


728 

4.10 

4.19 

2.84 

0.07 

l 3e 

4.04E-3 

0.20 

0.20412 


730 

4.13 

4.20 

2.84 

0.02 

l_4a 

1.85E-2 

0.05 

0.05153 


606 

1.83 

5.72 

4.31 

3.81 

l_4b 

1.22E-2 

0.08 

0.08276 


608 

4.49 

5.77 

4.36 

1.18 

l_4c 

9.85E-3 

0.10 

0.10341 


612 

5.06 

5.75 

4.34 

0.54 

l_4d 

6.65E-3 

0.15 

0.15512 


626 

5.53 

5.76 

4.35 

0.07 

l_4e 

5.00E-3 

0.20 

0.20665 


642 

5.60 

5.73 

4.33 

0.007 

TABLE III: Same as Table ED but for TOV model 2. 

Model 

Ao 

V 


0 

^•col 

Mah 

Madm 

Eq 

-fi’final 




1 

2 






2_la 

5.00E-5 

0.05 

0.04991 


750 

1.28 

1.31 

3.2E-5 

3.1E-5 

2_lb 

3.30E-5 

0.08 

0.07953 


750 

1.28 

1.31 

3.2E-5 

2.1E-5 

2_lc 

2.65E-5 

0.10 

0.09909 

0.09981 

750 

1.28 

1.31 

3.2E-5 

1.0E-5 

2_ld 

1.80E-5 

0.15 

0.14650 

0.14918 

750 

1.28 

1.31 

3.2E-5 

9.9E-6 

2_le 

1.35E-5 

0.20 

0.19802 

0.19927 

750 

1.28 

1.31 

3.2E-5 

2.7E-6 

2_2a 

6.30E-3 

0.05 

0.05009 


640 

1.28 

1.80 

0.50 

0.49 

2_2b 

4.12E-3 

0.08 

0.07971 

0.08025 

650 

1.54 

1.80 

0.50 

0.23 

2_2c 

3.33E-3 

0.10 

0.09859 

0.10008 

654 

1.64 

1.80 

0.50 

0.16 

2_2d 

2.25E-3 

0.15 

0.14936 

0.15009 

660 

1.67 

1.80 

0.50 

0.10 

2 2e 

1.70E-3 

0.20 

0.20017 


666 

1.76 

1.80 

0.50 

0.03 

2_3a 

1.50E-2 

0.05 

0.05117 


438 

1.36 

4.07 

2.84 

2.71 

2_3b 

9.82E-3 

0.08 

0.07918 

0.08170 

472 

3.20 

4.07 

2.83 

0.90 

2_3c 

7.95E-3 

0.10 

0.10197 


488 

3.39 

4.07 

2.83 

0.65 

2_3d 

5.36E-3 

0.15 

0.15314 


516 

3.98 

4.07 

2.83 

0.09 

2 3e 

4.04E-3 

0.20 

0.20412 


536 

4.01 

4.07 

2.83 

0.02 

2_4a 

1.85E-2 

0.05 

0.05170 


386 

1.58 

5.60 

4.30 

3.92 

2_4b 

1.22E-2 

0.08 

0.08276 


426 

4.32 

5.65 

4.35 

1.18 

2_4c 

9.85E-3 

0.10 

0.10341 


444 

5.05 

5.63 

4.33 

0.44 

2_4d 

6.65E-3 

0.15 

0.15512 


478 

5.40 

5.63 

4.33 

0.08 

2_4e 

5.00E-3 

0.20 

0.20678 


500 

5.48 

5.62 

4.32 

0.008 
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TABLE IV: Same as Table U but for TOV model 3. 


Model 

^0 



LO 

^col 

Mah 

Madm 

Eq 

-^final 




1 

2 






3_la 

5.00E-5 

0.05 

0.04991 


228 

1.09 

1.10 

3.2E-5 

3.1E-5 

3_lb 

3.30E-5 

0.08 

0.07971 


228 

1.09 

1.10 

3.2E-5 

2.6E-5 

3_lc 

2.65E-5 

0.10 

0.09928 

0.09982 

228 

1.09 

1.10 

3.2E-5 

1.4E-5 

3_ld 

1.80E-5 

0.15 

0.14758 

0.14936 

228 

1.09 

1.10 

3.2E-5 

1.1E-5 

3_le 

1.35E-5 

0.20 

0.19855 

0.19945 

228 

1.09 

1.10 

3.2E-5 

3.1E-6 

3_2a 

6.30E-3 

0.05 

0.05009 


226 

1.10 

1.60 

0.50 

0.49 

3_2b 

4.12E-3 

0.08 

0.07989 


228 

1.22 

1.60 

0.50 

0.36 

3_2c 

3.33E-3 

0.10 

0.09928 

0.10018 

228 

1.40 

1.60 

0.50 

0.16 

3_2d 

2.25E-3 

0.15 

0.14953 

0.15026 

228 

1.50 

1.60 

0.50 

0.09 

3 2e 

1.70E-3 

0.20 

0.19827 

0.20017 

228 

1.53 

1.60 

0.50 

0.04 

3_3a 

1.50E-2 

0.05 

0.05117 


220 

1.13 

3.91 

2.82 

2.71 

3_3b 

9.82E-3 

0.08 

0.07953 

0.08133 

226 

2.94 

3.91 

2.82 

0.91 

3_3c 

7.95E-3 

0.10 

0.10197 


226 

3.22 

3.91 

2.82 

0.61 

3_3d 

5.36E-3 

0.15 

0.15314 


228 

3.72 

3.90 

2.81 

0.11 

3 3e 

4.04E-3 

0.20 

0.20412 


228 

3.82 

3.91 

2.82 

0.03 

3_4a 

1.85E-2 

0.05 

0.05170 


218 

1.23 

5.35 

4.28 

4.02 

3_4b 

1.22E-2 

0.08 

0.08276 


226 

4.15 

5.39 

4.32 

1.15 

3_4c 

9.85E-3 

0.10 

0.10341 


226 

4.82 

5.35 

4.30 

0.43 

3_4d 

6.65E-3 

0.15 

0.15512 


228 

5.17 

5.38 

4.31 

0.08 

3_4e 

5.00E-3 

0.20 

0.20682 


228 

5.29 

5.37 

4.29 

0.01 


TABLE V: Same as Table ED but for TOV model 4. 

Model 

Ao 


to 

1 2 

^-col 

Mah 

Madm 

E 0 

E final 

4_la 

5.00E-5 

0.05 

0.04991 


72 

0.75 

0.80 

3.2E-5 

3.1E-5 

4_lb 

3.30E-5 

0.08 

0.07989 


72 

0.75 

0.80 

3.2E-5 

3.1E-5 

4_lc 

2.25E-5 

0.10 

0.09963 


72 

0.75 

0.80 

2.3E-5 

1.9E-5 

4_ld 

1.80E-5 

0.15 

0.14882 

0.14973 

72 

0.75 

0.80 

3.2E-5 

1.4E-5 

4 le 

1.35E-5 

0.20 

0.19927 

0.19963 

72 

0.75 

0.80 

3.2E-5 

7.5E-5 

4_2a 

6.30E-3 

0.05 

0.05009 


72 

0.75 

1.29 

0.50 

0.49 

4_2b 

4.12E-3 

0.08 

0.08007 


72 

0.77 

1.29 

0.50 

0.47 

4_2c 

3.33E-3 

0.10 

0.09982 

0.10035 

72 

0.89 

1.28 

0.49 

0.35 

4_2d 

2.25E-3 

0.15 

0.14829 

0.15008 

72 

1.09 

1.29 

0.50 

0.19 

4_2e 

1.70E-3 

0.20 

0.19908 

0.20035 

72 

1.24 

1.29 

0.50 

0.04 

4_3a 

1.50E-2 

0.05 

0.05117 


72 

0.84 

3.60 

2.80 

2.72 

4_3b 

9.82E-3 

0.08 

0.08078 

0.08168 

72 

2.61 

3.59 

2.79 

0.91 

4_3c 

7.95E-3 

0.10 

0.10214 


72 

2.63 

3.60 

2.80 

0.88 

4_3d 

5.36E-3 

0.15 

0.15314 


72 

3.38 

3.59 

2.79 

0.11 

4 3e 

4.04E-3 

0.20 

0.20412 


72 

3.48 

3.60 

2.80 

0.04 

4_4a 

1.85E-2 

0.05 

0.05182 


72 

0.86 

5.03 

4.24 

4.08 

4_4b 

1.22E-2 

0.08 

0.08029 

0.08276 

72 

3.53 

5.04 

4.29 

1.45 

4_4c 

9.85E-3 

0.10 

0.10341 


72 

4.30 

5.06 

4.27 

0.64 

4_4d 

6.65E-3 

0.15 

0.15512 


72 

4.82 

5.06 

4.27 

0.10 

4_4e 

5.00E-3 

0.20 

0.20681 


72 

4.93 

5.05 

4.26 

0.01 
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FIG. 4: Left column: Time evolution of the scalar field with initial energy Eo = 3.2 x 10~ 5 for the TOV model 2 and three 
different scalar field masses g, = {0.05,0.1,0.2} (from top to bottom). The inset shows a magnified view of the initial 2000 
units of time in the evolution. Right column: Fourier transform of the evolution of the scalar field, to refers to the real part 
of the frequencies. The units in the vertical axis are arbitrary and are normalized by the amplitude of the fundamental mode. 
Each curve corresponds to a different extraction radius, as indicated in the legend of the top-left plot. 
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FIG. 5: Time evolution of the scalar field with initial energy 
Eq = 4.3 for the TOV model 2 and three different scalar field 
masses /r = {0.05,0.1,0.2} (from top to bottom). The inset 
shows a magnified view of the initial 2000 units of time in the 
evolution. 


FIG. 6: Time evolution of the scalar field with initial energy 
Eo = 4.3 for the TOV model 4 and three different scalar field 
masses /r = {0.05,0.1,0.2} (from top to bottom). The inset 
shows a magnified view of the initial 2000 units of time in the 
evolution. 
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FIG. 7: Time evolution of the maximum (central) rest-mass 
density of the TOV model 1 for different initial values of the 
amplitude of the scalar field pulse with mass fj, = 0.05. Notice 
that the black and red solid lines overlap. 


alyze the results of the simulations by extracting a time 
series for the scalar field amplitude at a set of observation 
points located at fixed radii r ex t (typically at three radii 
located at r ex t = 100,200 and 300). To identify the fre¬ 
quencies at which the scalar field oscillates we perform a 
Fast Fourier transform after a given number of time steps 
and obtain the corresponding power spectrum. The very 
existence of these frequencies for the long-term simula¬ 
tions we perform is the main indicator we use to demon¬ 
strate the presence of quasi-bound states. 

The main results of our simulations are summarized 
in Tables HIHVl These four tables contain for each TOV 
model, the different choices for the scalar field mass /q 
the initial amplitude of the pulse Ao, the Arnowitt-Deser- 
Misner (ADM) mass Madm of the whole system, defined 
as 


lim (f djipdS-i 
r ^°°Js 

and the initial energy of the scalar field E 0 


Madm = — x — — 

Z7T r—¥oo / g 


E = 


/•OO 

/ ‘ 

J2M 


£ SF dV 


(39) 


(40) 


where £ SF is defined in Eq. (181. As the second-to-last 


column of each table shows, we vary the initial ampli¬ 
tude for the different /i in order to keep the initial scalar 
field energy constant for four different scenarios, from the 
test-field approximation to the non-linear regime where 
the scalar field is self-gravitating. This is presented as the 
four subdivisions of each table. The oscillations frequen¬ 
cies of the scalar field are reported in columns 4 and 5 for 
the fundamental mode and first overtone, respectively. 

Figures [4][6] show the time evolution of the scalar 
field amplitude extracted at three different radii, r ex t = 


{100,200,300}. In the left panels of Fig. [4] we plot the 
evolution for the initial scalar field energy Eq = 3.2 x 10~ 5 
for the TOV model 2 and for three scalar field masses, 
namely /r = {0.05,0.1,0.2} from top to bottom. The 
three panels show a moderately rapid decay of the am¬ 
plitude of the oscillations of the scalar field. However, the 
field amplitudes by the end of the simulations have not 
totally vanished, especially for the fj, = 0.05 model shown 
in the top panel. Correspondingly, the right panels of 
Fig.[4]present the power spectra, showing the well-defined 
oscillation frequencies of the quasi-stationary states. It 
is worth stressing that the values of the frequencies do 
not depend on the location of the observer. 

Similar time evolution plots are displayed in figure [5] 
but, in this case, for the larger initial scalar field en¬ 
ergy of our sample, E 0 = 4.3 (keeping the same values of 
H and TOV model as in Fig. [4]). The evolution for the 
same initial parameters corresponding to the TOV model 
4 can be found in figure [6] The comparison among these 
time evolution plots shows that for the smallest scalar 
field mass, jj = 0.05, the differences between the evo¬ 
lutions for the different TOV models are more evident. 
The TOV model 2 shows an increase in the amplitude of 
the oscillations at the end of the simulation due to the 
presence of two close frequencies giving rise to a beating 
pattern, while TOV model 4 does not show this behavior. 
However, for larger scalar field masses, the differences in 
the time evolutions become qualitatively less important. 

The power spectra obtained from the Fourier trans¬ 
forms of these time series show a set of distinct frequen¬ 
cies, exemplified in the right panels of Fig.[4]for a partic¬ 
ular set of models. The frequency resolution is inversely 
proportional to the simulation time and hence the longer 
the runs the more accurate the frequencies. The presence 
of these frequencies confirms the existence of long-lived 
spherical quasi-bound states of scalar fields surrounding 
black holes even if these are formed dynamically through 
the gravitational collapse of spherical stars. During the 
non-linear evolution of the system, part of the scalar field 
survives the gravitational collapse and remains in the 
form of trapped states. These results extend the valid¬ 
ity of our previous model, put forward in [32] and Paper 
1 , in the linear and non-linear regimes respectively, to a 
significantly more dynamical scenario. 

In columns 4 and 5 of Tables |HffV| we show the frequen¬ 
cies u> of the scalar field oscillations. These frequencies 
refer only to the real part of the dominant modes and are 
labelled according to their strength on the power spec¬ 
trum. For some models we are able to obtain several fre¬ 
quencies of oscillation. Tables m also report the mass 
of the apparent horizon (AH) at the end of the simula¬ 
tion Mah, defined as Mah = -\/A/ 16 tt, where A is the 
area of the AH, the time of collapse at which the AH is 
first found, f co i, and the scalar field energy Eft na \ at the 
end of the simulation. In figure [7] we plot the time evo¬ 
lution of the central rest-mass density of the TOV model 
1 for different values of the scalar field amplitude. This 
figure shows that the time at which the TOV collapses 
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decreases monotonically when the initial scalar field am¬ 
plitude increases. We note that t co \ for the case of the 
test-field scalar field energy (A 0 = 5 x 10 -5 ) coincides 
perfectly with the corresponding time for the case with 
no scalar field affecting the density evolution of the TOV 
star (as both lines overlap in the figure). 


From Tables [Iiyv| we see that the time of collapse shows 
a clear dependence with the initial energy of the scalar 
field E 0 only for the less compact TOV models of our 
sample, models 1 and 2 for which the values of the com¬ 
pactness parameter M/R are 0.008 and 0.015, respec¬ 
tively. This time becomes shorter the larger Eg. Model 
3, with M/R = 0.025, has significantly shorter values of 
f co i and, moreover, t co i barely shows any dependence on 
Eq. Finally, for the TOV model 4 with M/R = 0.042 
the time of collapse is the shortest of the whole sample 
and constant, t co \ = 72 irrespective of Eq. This result is 
explained by the fact that in this case the star collapses 
to form a black hole even before the scalar field pulse can 
reach the collapsing star. Therefore, the presence of the 
scalar field does not influence the collapse dynamics at 
all, even for the highest scalar field energies we consider 
in our sample. 

The oscillation frequencies for the same scalar field pa¬ 
rameters vary slightly from one TOV model to another 
due to the different initial mass of the black hole when 
it is formed. For most cases we find that the smaller the 
black hole mass the higher the frequencies. This trend, 
however, does not apply for the larger energy models, for 
which the frequencies are essentially identical for all the 
TOV stars. This behaviour might be simply a numerical 
artifact due to the lack of sufficient time resolution, as 
with increasing evolution time the frequencies could be 
better resolved and differences may appear. 

In general our results show that when we consider large 
enough values of the scalar field mass /q a significantly 
large amount of the initial scalar field energy falls onto 
the black hole, and this amount increases with p, (see 
last columns of Tables 00. This behaviour is found for 
all initial scalar field energies considered. On the other 
hand, the amount of scalar field absorbed by the black 
hole changes depending on the TOV model, that is, on 
the initial mass of the black hole when it is formed. We 
find that it is slightly smaller for the models for which 
the TOV mass is smaller (and the compactness of the 
star is correspondingly larger). For the models with the 
largest mass ^ and largest amplitude Aq , the energies at 
the end of the simulations are just a tiny fraction of the 
initial values (see last row of last columns of Tables n]|V I. 
Even if long-lived quasi-stationary states can be found in 
such cases, their energies may be too small to question 
the actual survival of the scalar field. Nevertheless, for 
smaller values of /r this trend is reversed and only a small 
fraction of the initial energy of the scalar field is absorbed 
by the black hole. 


Our results cannot be directly extrapolated to the case 
of a typical 1O 6 M 0 SMS collapsing to a black hole. How¬ 
ever, we may still put forward an educated guess of what 



0 10000 20000 30000 

Time 


FIG. 8: Time evolutions of the scalar field energy from the 
initial value Eq = 3.2 x 10~ 5 (top panel ) and from the ini¬ 
tial value Eo ~ 4.3 (bottom panel). The evolutions are 
shown for the TOV model 2 and for all scalar field masses 
p = {0.05,0.08,0.1,0.15,0.2}. The insets show magnified 
views around the time of collapse in the corresponding evo¬ 
lutions. The times at which the black holes are formed are 
displayed as vertical dashed lines. The different colors are 
related to the different scalar field masses as indicated in the 
legend of the top panel. 


actually might happen. On the one hand, SMS are less 
compact ( M/R ~ 10 -3 ) than our TOV models and their 
gravity is almost Newtonian. Taking this into account, 
if the quasistationary states survive to the collapse of a 
more compact star, it could be expected that they will 
also persist after the collapse of a SMS. On the other 
hand, SMS collapse occurs at f col ~ 10 5 s [5], which 
in our units is about t co \/M ~ 1000 — 40000, depend¬ 
ing on the compactness of the star. Our results yield 
tcoi/M ~ 90 — 1750, and the simulations show that the 
scalar field survives at least for two to three orders of 
magnitude beyond the collapse timescale. In the context 
of a SMS, quasistationary states may then exist for long 
times. 

The time evolution of the scalar field energy corre- 
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Time 

FIG. 9: Time evolution of the mass of the apparent horizon 
A4 ah and of the scalar field energy for the initial Eo ~ 4.3 
and the scalar field masses p = {0.05,0.1,0.2} from top to 
bottom corresponding to the TOV model 2. The inset shows 
a magnified view of initial 2000 in the evolution. The time 
at which the black hole is formed is displayed as a vertical 
dashed line. 


sponding to the TOV model 2 and for all scalar field 
masses is shown in figure [8] The top panel corresponds 
to Eo = 3.2 x 10 -5 and the bottom panel to E 0 ~ 4.3 (the 
precise values are given in the second to last column of 
Table |Hl} models 2.4a-2.4e). In the case of the test-field 
regime, the energy decay is smooth and varies monotoni- 






p 


FIG. 10: Top panel: Frequency u as a function of the scalar 
field mass p for the models shown in Table |III| The semi- 
analytic relation given by Eq. (41 1 is also plotted for M = 
{1.31,1.5, 2.0} (solid lines). The semi-analytic values for M = 
1.31 corresponds to the cyan line. Bottom panel, same as the 
top panel but for the models shown in Table [V] The semi- 
analytic relation is plotted for M = {0.8,1.5, 2.0}. The semi- 
analytic values for M = 0.80 corresponds to the cyan line. 


cally with the increasing mass. In the non-linear regime, 
there is a significant leap (better seen in the inset of the 
bottom panel) in the amount of scalar field energy due 
to the initial accretion onto the black hole when it is 
formed. The decay of the energy after black hole forma¬ 
tion (signalled by the dashed vertical lines in the figure) is 
in general fairly abrupt for all scalar field masses except 
for p = 0.05. After the scalar field first interacts with 
the black hole, the energy reaches a phase during which 
the accretion rate is importantly reduced, allowing to the 
survival of the scalar field for long timescales (but again 
notice the drastic reduction in the scalar field energy for 
E 0 ~ 4.3 and p = 0.15 and p = 0.2). This response of 
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the system happens faster the larger the initial energy of 
the scalar field. 

In Fig. [9]we plot the evolution of the scalar field energy 
Esf (red curves) and the mass of the apparent horizon 
Mah (black curves) for the TOV model 2 and the largest 
initial energy (Eq ~ 4.3) for three scalar field masses, 
H = {0.05,0.1,0.2} from top to bottom. The scalar field 
energy remains mostly constant until the time of collapse, 
when it begins to decrease (see the corresponding insets). 
The decrease is more significant and faster the larger the 
scalar field mass. Similarly, the mass of the black hole 
starts to increase once the apparent horizon (discontinu- 
ally) appears, as a response to the accretion of the scalar 
field. It becomes clear from this figure that the amount 
of the scalar field lost in the process is absorbed by the 
black hole, becoming the only contribution to its growth 
when the fluid mass has been completely accreted by the 
black hole. The evolutions of both quantities are almost 
perfect mirror images from one another. 

By using a matching technique, Furuhashi and Narnbu 
showed analytically in (33j that in the limit 1 the 

real part of the frequency of quasi-bound states, w q b s , 
depends on the mass parameter /i as 


Re(wqb s ) ~ /x 


l--(Mp) 2 


(41) 


In Fig. [To] we show a plot of the frequencies of all the con¬ 
figurations of Table III and [v] as a function of /x (similar 
behaviour is found for the other models). Our numerical 
results are plotted with non-solid lines. We also consider 
in this plot the adiabatic approximation for the frequency 
given by Eq. (411 for different values of the mass of the 


black hole M. These analytic results are plotted with 
solid lines in the figure. 

For the TOV model 4 (bottom panel), we find that for 
the smallest value of Aq, indicated in the legend in the 
top panel of the figure, the corresponding semi-analytic 
curve (black dotted line) is indistinguishable from the one 
of the test field approximation for M = 0.8 (cyan solid 
line) showing consistency with the test-field approxima¬ 
tion and the limit M/i< 1. For the TOV model 2 (top 
panel), while the results are similar, for the larger scalar 
field mass the adiabatic approximation is not matched 
even in the test-field approximation for M = 1.31 (cyan 
solid line). For all the models, we find that for greater 
values of Aq, the results from the non-linear approach 
not only deviate from the adiabatic approximation but 
they also follow the opposite trend: the frequency is 
greater than the frequency of the test-field limit. This 
discrepancy may be due to the violation of the condition 
M/i <C 1 used in the analytical approximation and due 
to a breakdown of the adiabatic approximation. 

These results agree with our findings in Paper 1. When 
the initial amplitude is increased and the model is far 
from the test-field regime, the frequencies become higher. 
A similar trend is found when the scalar field mass is 
small, the frequencies becoming closer to those of the 
test-field regime. However, in this case we find that for 


the larger amplitudes the frequencies do not differ as 
much as in Paper 1 from the frequencies matching the 
adiabatic approximation. This may be due to the fact 
that depending on the initial profile and position of the 
scalar field pulse, there can be a variable amount of scalar 
field that ends up being part of the quasi-bound states. 


C. Comparison with a BH-scalar field system 


In order to further understand the influence of the col¬ 
lapse of the stars in the evolution of the scalar field, we 
also evolve an additional set of 12 models corresponding 
to purely Schwarzschild-like black hole spacetimes. These 
black holes have initial masses equal to those of the dif¬ 
ferent TOV stars considered in the simulations we have 
just discussed. Likewise, the initial distributions of the 
scalar field clouds are the same as in the previous sim¬ 
ulations. Therefore, the scalar field interacts from the 
beginning with an already existing black hole. This is 
the same situation we analyzed in Paper 1, only the ini¬ 
tial black hole masses have now been adapted to match 
those of the TOV stars. We choose only three scalar field 
masses, /x = {0.05,0.1,0.2}, and the largest value of the 
initial energy, Eq ~ 4.3. The grid resolution, time-step, 
and final evolution time remain unchanged. Table |VI| 
summarizes the results corresponding to this scenario. 

In this set of simulations, as we showed in Paper 1, part 
of the scalar field is immediately accreted onto the black 
hole while part forms the quasi-bounds states. The de¬ 
termination of the oscillation frequencies of those states 
reveals that, for /i = 0.1 and 0.2, some of the frequen¬ 
cies are equal to those of the collapsing star scenario. 
On the other hand, the final scalar field energy is, for 
the majority of models, smaller than that for the TOV 
setup but the difference is reduced either if the scalar 
field mass is increased or if the initial black hole mass 
is lowered. This discrepancy, most noticeable for TOV 
model 2, can be due to the fact that for the TOV a larger 
fraction of the scalar field may escape before the black 
hole forms and starts accreting. There is almost no differ¬ 
ence between the two setups for TOV model 4 because, 
as mentioned before, the time of collapse is short enough 
to form a black hole before most of scalar field reaches 
the star. We also note that for the black hole models 1 
and 2 and /x = 0.05 the amount of scalar field that falls 
onto the black hole is larger than for the other models 
and scalar field masses. As a result, the mass of the fi¬ 
nal black holes increases significantly and the oscillation 
frequencies change. The dynamics of the scalar field in 
the BH and TOV setups are markedly different for this 
particular set of models. 

This is clearly visible in Fig. [TT] where we plot the evo¬ 
lution of the scalar field energy for the initial parameters 
described in Table [Vj for both the TOV and BH models. 
As the mass of the scalar field increases from [i = 0.05 
(left column of panels) to /x = 0.2 (right column of pan¬ 
els) the evolution of the energy becomes progressively 
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TABLE VI: Same as Tables |U|V| but for evolutions of scalar fields in spacetimes containing initially a Schwarzschild black hole. 
Only the fundamental frequency is shown. 


Model 

A 0 

A 1 

UJ 

(Mah)o 

(TIah) final 

A/adm 

E 0 

E final 

l_4a_bh 

1.85E-2 

0.05 

0.05116 

1.43 

3.70 

5.73 

4.31 

1.88 

l_4c_bli 

9.85E-3 

0.10 

0.10341 

1.43 

5.32 

5.75 

4.34 

0.31 

l 4e bli 

5.00E-3 

0.20 

0.20681 

1.43 

5.62 

5.75 

4.33 

0.007 

2_4a_bh 

1.85E-2 

0.05 

0.05170 

1.31 

2.67 

5.59 

4.30 

2.88 

2_4c_bli 

9.85E-3 

0.10 

0.10341 

1.31 

5.18 

5.62 

4.33 

0.34 

2 4e bh 

5.00E-3 

0.20 

0.20681 

1.31 

5.44 

5.61 

4.32 

0.008 

3_4a_bh 

1.85E-2 

0.05 

0.05170 

1.10 

1.34 

5.36 

4.28 

3.93 

3_4c_bh 

9.85E-3 

0.10 

0.10341 

1.10 

4.87 

5.39 

4.30 

0.41 

3 4e bli 

5.00E-3 

0.20 

0.20681 

1.10 

5.26 

5.38 

4.29 

0.01 

4_4a_bh 

1.85E-2 

0.05 

0.05188 

0.80 

0.86 

5.04 

4.24 

4.08 

4_4c_bli 

9.85E-3 

0.10 

0.10341 

0.80 

4.32 

5.06 

4.27 

0.64 

4_4e_bh 

5.00E-3 

0.20 

0.20681 

0.80 

4.93 

5.05 

4.26 

0.01 


similar, and the presence of a black hole from the start 
or at later times is unimportant for the dynamics. The 
same trend is found as the initial black hole mass or TOV 
mass is reduced (see bottom row of panels in Fig. 111. In 
this case the dynamics also becomes identical irrespec¬ 
tive of the scalar field mass. The major differences are 
found in the upper left corner of the composite of panels 
in Fig. 0i-e- when [i is the smallest and the black hole 
and TOV initial masses are largest. In any event, the 
influence of the gravitational collapse on the dynamics 
of the scalar field and on the eventual presence of quasi¬ 
bound states is not as potentially threatening as we might 
have expected before conducting this research. The ini¬ 
tial (rapid) phase of accretion by an already formed black 
hole seems to overcome any dynamical effect from the col¬ 
lapse of the star which might have affected the survival 
of the scalar field, even for the smaller scalar field masses 
and the larger initial black hole masses considered in our 
study, where different evolutions are observed. 


V. CONCLUSIONS 


In this paper we have discussed the feasibility of the ex¬ 
istence of quasi-stationary states of self-gravitating scalar 
fields in the form of clouds (or hairy “wigs”) around col¬ 
lapsing, black-hole-forming, stars. For our analysis we 
have solved numerically the Einstein-Klein- Gordon cou¬ 
pled system of equations further coupled to the Euler 
equations of general relativistic hydrodynamics. Our nu¬ 
merical approach has been limited to spherical symmetry 
and has been based on spherical coordinates and a PIRK 
numerical scheme for the time update. This research 
extends our previous work where the same issue was in¬ 
vestigated for black holes already existing at t = 0 to 
a much more dynamical scenario where the black holes 
are formed as the result of the gravitational collapse of 
unstable (polytropic) stars. 


In our non-linear study we have performed a large 
number of (second-order) accurate and long-term stable 
simulations of dynamical non-rotating collapsing stars, 
with an EOS corresponding to SMSs, surrounded by self- 
gravitating scalar fields. The models of our sample have 
been suitably parameterized to span a broad range of 
scalar field cases, from the test-field regime to the fully 
non-linear regime, along with different initial masses for 
the star. We have found that in the case of highly dy¬ 
namical spacetimes there are states which closely resem¬ 
ble the quasi-bound states in the test-field approximation 
of [2D] [32], surviving the stellar collapse and black hole 
formation. 

By performing a Fourier transform of the time series 
of our numerical data we have been able to characterize 
the scalar field states by their distinctive oscillation fre¬ 
quencies. It has been found that the scalar field oscillates 
with different well-defined frequencies, and that its time 
evolution produces in most of our models a character¬ 
istic beating pattern due to the non-linear combination 
of two oscillations with close frequencies. Such scalar 
field oscillations may have important imprints in a num¬ 
ber of astrophysical scenarios such as gravitational wave 
astronomy. In particular m inn have shown recently 
that the frequency of oscillation of the scalar field leaves 
a distinctive signature in the late time behaviour of the 
gravitational wave signal. 

We have also compared the evolutions of the scalar field 
in collapsing, black-hole-forming spacetimes with space- 
times that already contain a non-rotating black hole. 
This comparison has provided a better understanding of 
the influence of the stellar collapse on the scalar field 
dynamics. Our results have shown that the differences 
are mainly due to the initial (or lack of) accretion phase 
onto the black hole that contributes to the growth of its 
mass and to the different dynamics, in the case of self- 
gravitating scalar fields. 

Furthermore, our results indicate that for certain val- 
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FIG. 11: Comparison of the time evolution of the scalar field energy for the collapsing TOV star models and for Schwarzschild 
black holes with the same initial mass. The initial energy density is Eq ~ 4.3. The insets show magnified views in logarithm 
scale of the initial 10000 units of time in the evolution. Left column: [i = 0.05. Middle column: n = 0.1. Right column: fi = 0.2. 
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ues of /z, the formation of the black hole during the col¬ 
lapse of the star enhances the chances of the scalar field 
to remain in a quasi-bound state. This result may have 
important implications for models of black hole forma¬ 
tion together with dark matter halos in the early epochs 
of the Universe. For instance, large galaxies in the scalar 
field model can already be formed at high redshift 
The collapse of the ~ IMq stars considered in this 
work and the subsequent formation of the black holes 
modify the dynamics of the scalar field, but for the range 
of TOV masses this type of events has very short dura¬ 
tion. While our results cannot be directly extrapolated 
to the case of a typical 1O 6 M 0 SMS collapsing to a black 
hole, it could perhaps be expected, on the basis of their 
smaller compactness and almost Newtonian gravity, that 
quasistationary states might still persist after the collapse 
of a SMS star. Actual numerical simulations are however 
necessary to validate the viability of this scenario and 
the survival of quasi-bound states of scalar fields around 
SMBH for cosmological timescales. In addition to such 
simulations it is worth to study the evolution of scalar 
fields around accreting black holes during significantly 


longer evolutionary times than those attempted in im. 
a task we defer for the near future. 

Finally, we note that the dynamical scenario presented 
here can also be generalized to rotating stars where the 
final state may be a rotating Kerr black hole surrounded 
by a scalar field cloud yielding a hairy black hole solu¬ 
tion mm- It has been argued that these solutions may 
act as dynamical attractors in a highly dynamical pro¬ 
cess |38| and thus simulations as the ones performed here 
would be able to describe their formation. 
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